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Dynamical symmetry breaking in three-dimensional QED with N fcrmion flavours is consid- 
ered at finite temperature, in the large N approximation. Using an approximate treatment 
of the Schwinger-Dyson equation for the fermion self-energy, we find that chiral symmetry 
is restored above a certain critical temperature which depends itself on N. We find that 
the ratio of the zero-momentum zero-temperature fermion mass to the critical temperature 
has a large value compared with four- fermion theories, as had been suggested in a previous 
work with a momentum-independent self-energy. Evidence of a temperature-dependent crit- 
ical N is shown to appear in this approximation. The phase diagram for spontaneous mass 
generation in the theory is presented in T — N space. 



1 Introduction 



Quantum Electrodynamics in (2+1) dimensions (QED3) has attracted considerable interest in 
the last few years [Q] - |5| . One reason for this is that it provides a simple setting for the study of 
dynamical chiral symmetry breaking which is important for theories such as QCD. Furthermore 
QED3 also appears to be relevant to some long-wavelength models of certain 2D condensed 
matter systems, including high-T c superconductors [Q, 0]. At zero temperature a number of 
studies have shown that chiral symmetry is dynamically broken in QED3. Using the leading 
order in the 1/N expansion of the Schwinger-Dyson (SD) equations Appelquist et al. Q showed 
that the theory exhibited a critical behaviour as the number iV of fermion flavours approached 
N c = 32 /it 2 ; that is, a fermion mass was dynamically generated only for N < N c . Qualitatively 
the same behaviour was found by Nash who included 0(1/N 2 ) corrections. As against this, 
Pennington and collaborators ||, adopting a more general non-perturbative approach to the 
SD equations, found that the dynamically generated fermion mass decreased exponentially with 
N, vanishing only as N — » 00 (as originally found by Pisarski Q using a simplified form of of 
the SD equation). On the other hand, an alternative non-perturbative study by Atkinson et al. 
|| suggested that chiral symmetry is indeed unbroken at sufficiently large N. We should also 
mention that Pisarski |5j has used the renormalization group approach and the e expansion to 
argue that chiral symmetry remains broken for all N, but this result has to be interpreted with 
some caution as the relevant value of e is 1. The theory has also been simulated on the lattice 
Q and the results appear to be consistent with the existence of a critical N as predicted in 
the analysis of ref Q. On the other hand, because Monte-Carlo simulations are performed on 
lattices of finite size L and typically cannot detect mass scales less than the IR cutoff scale 1/L, 
the persistence of an exponentially small fermion mass for large N, as suggested in Q, cannot 
be ruled out. 

The extension of the above type of analysis to finite temperature is extremely important, 
and highly relevant to either application already mentioned. The most obvious question is 
whether there is a critical temperature T c above which chiral symmetry is restored. This was first 
answered affirmatively by Kocic J10|| , using a very simple approximation to the finite temperature 
SD equations, in which the entropy of the fermions was not fully taken into account. An 
improved calculation was made by Dorey and Mavromatos []1 1[] , based on the finite temperature 
Schwinger-Dyson equations, to leading order in 1/N. These authors found that the ratio r of 
twice the zero-temperature mass to the critical temperature was approximately independent of 
N and much larger (r ~ 10) than the value obtained in (BCS-like) four fermion theories (r at 
most 3.5). The latter result could be relevant to an understanding of this ratio in the high-T c 
superconductors. In ref pT| , however, the (leading order in 1/N) SD equation for the fermion 
self-energy function £(p) was considerably simplified by making the assumption ||] that X was 
in fact a constant, independent of p. It is clearly necessary to be assured that the results of [11 
are essentially independent of this assumption. 

In the present paper we extend and complete the analysis of Ref. [11|. We again start 
from the Schwinger-Dyson equation for X(p), to leading order in 1/N, but keep the full momen- 
tum dependence of X. Following the approach of [fLl] ], and of similar studies of chiral symmetry 
restoration in four dimensions we further simplify the SD equation by adopting an instan- 
taneous approximation for the kernel. This corresponds to retaining only the part of the kernel 
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corresponding to the static interaction between charges. We find that the main conclusions of 



Ref. [11] are indeed robust -in particular, the ratio r has almost the same value. 

A second question concerns the existence of a critical number of flavours N c , above which 
chiral symmetry would be restored. Within the 1/N approach, there is no N c if £ is approx- 
imated by a constant [pi, as already noted, a result which naturally persists at finite T in the 



approximation |ll[ ; but with ap-dependent £ we find a (temperature dependent) N c , analogous 
to Ref. ||]. We are thus able to present the resulting phase structure of the model, as a func- 
tion of the variables N and T: there is a single critical line separating the region of low (N, T) 
where S / from the region of high (N,T) where £ = (see Figure 7 below). As indicated 
above, we are aware that criticism have been raised || concerning the reliability of the 1/N 
approach for calculating fermion masses £ well below the intrinsic scale set by the dimensionful 
coupling a (as happens in this case) and it may be that in the exact theory no sharply defined 
N c exists (though it would seem very difficult to establish this conclusively, given that some 
approximation to the complete set of SD equations has to be made). However, at the very 
least, the lattice results at T = strongly indicate the existence of an effective critical region, 
in which the dynamical fermion mass undergoes a rapid crossover from a regime of large values 
to one of very small values as N is increased. We shall take the view here that our 1/iV-based 
phase diagram should give a qualitatively correct picture in this sense -and that in any case 
there is merit in obtaining as complete an analysis as possible of one definite model at finite 
temperature. It would clearly be of interest to investigate the effect, at finite temperature, of 
the non-perturbative vertex structure advocated in Ref. Q. 



2 Momentum-dependent self-energy (gap) equation at finite tem- 
perature 

The Lagrangian for massless QED3 with N flavours is 

C = -\UuF u + Mi^-e^ (1) 

where is the vector potential, i = 1,2, ...N, and a reducible four dimensional Dirac algebra 
has been chosen so that ([!]) has a continuous chiral symmetry. The conventions of [|| will be 
adopted throughout. 

The Schwinger-Dyson equation for the fermion propagator at non-zero temperature &bT = 
/3 _1 is given by 

+00 „ j2k. 

S F 1 ( P0 ,P,P) = SP- I (p)-^ £ —-^rS F (k ,K,P)A llu (qo,Q,P)r'; (2) 

' n=— 00 ' 

P = (Po,p) P = IpI Po = (2m + l)7r//3 

k = (k ,k) K=\k\ k = (2n + l)vr//3 (3) 

<?=(<7o,q) Q = |q| = |p-k| q = 2(m - n)vr//? 



where 
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As stated above, we truncate (|2|) by working at leading order in 1/N, in which case T u is 
replaced by its bare value e^ v and A^ u by the 0(1/N) propagator shown in Figure 1, in which 
the fermions are massless (the massless vacuum polarisation loop already softens the photon 
propagator jjj , and the exact form of this softening does not qualitatively change the behaviour 
of the fermion propagator ||). We shall work in Landau gauge and assume that the wave 
function renormalisation can be neglected to leading order in 1/N, so that Sp 1 = $ + S m (P, (3) 
(note however the criticism made of this step by Pennington and Walsh Q). The trace of 
equation (^) then yields a closed integral equation for T, m : 

with 

A(<?o,Q,/?) = -Tr^A^Q^hl (5) 

and a = Ne 2 . We now follow Ref. [11] in retaining only the fj, = v = component of A^ u at 
zero frequency and thus set, 

vfo.q.fl = 7gjg%i * ( 6 ) 

Q z + U (Q,/3) 

where [11| 

n (Q,/?) = — R tdx ln(2 cosM^J^l-x))). (7) 
np Jo 2 v 

In this approximation S m (P) becomes frequency independent and the summation over n 
in (U) can be performed analytically yielding 

e( p,„) _ - "i^y, ( 8 ) 



siv^y Q 2 + n (g,/3) v^Ts^cO) 

The main purpose of this paper is to present the detailed analysis of this temperature- and 
momentum- dependent gap equation. 



3 Numerical procedure and results 

Equation (||) involves a two-dimensional integral over k, while Ho in (0) involves Q = |p— k| inside 
a further integral. The analogous T = equation has frequently been simplified |l^] by replacing 
Q by max(K, P), rendering the angular integral trivial. We do not make this approximation 
here, but we have found it very convenient to adopt an excellent analytic approximation to (0) 
(correct to about 1.5% at worst) which is provided by the expression 

n (Q,/5) = ^ 



16 In 2 

Qp H exp ' 



7T 



7T 



16 In 2 



QP) 



(9) 
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which incorporates the correct limiting behaviour as either Q or [3 tends to zero or infinity. In 
particular, as Q — > 0, 

n (Q,/?)^^ (io) 

which exhibits the thermal screening noted before Jll[| . 

As regards to the integral over if = |k| in (||), this of course extends to K — > oo in 
principle. For numerical purposes the upper limit would normally be replaced by some cutoff 
parameter A, chosen to be sufficiently large that further increase of it makes no difference. In 
the present case however, the dimensionful parameter a provides a natural scale: in particular, 
Appelquist et al. Q noted that the integral in the corresponding T = equation was rapidly 
damped for momenta greater than a, so that effectively A ~ a. We have found that the same 
is true for equation (||) at finite temperature, and thus from now on we shall work with (Js[) cut 
off at A = a, and present our results in terms of the scaled quantity (E/a), momenta being also 
scaled by a. 

We have used a numerical algorithm to solve equation <^j. An iterative procedure has 
been followed to find a solution for £(P, (3). As usual the algorithm only converged if the input 
function was sufficiently close to the true solution. The iteration was started by adopting the 
approximation Q ~ max(K, P) and solving the resulting version of (|8|) for a large value of (3. 
This solution then provided the input trial function for the true equation (||) at the same low 
temperature. When the input and the output of the integral equation agreed within a 2% of 
difference the iterative sequence was stopped, and the output was considered as the solution to 
equation (^) for the defined temperature. f3 was then incremented downwards in small steps 
using the solution of (||) for the previous (3 as the input function for the next. In this way we 
were able to obtain the dependence of £ on (3, as well as on p. 

The scaled dynamical mass as a function of scaled momentum is shown in Figure 2 for 
A = 1, at various temperatures. We notice that the mass remains constant for a wide range of 
momenta up to roughly P ~ S(P = 0) and then rapidly drops to zero. We also notice that as 
expected the mass decreases with rising temperature. The rate at which this occurs grows as 
we approach the critical temperature T c , above which the mass vanishes (and which will depend 
on TV). This is illustrated in Figure 3, which shows the zero momentum mass S(0, (3) /a as a 
function of the scaled temperature ksT/a for A = 1,1.5. The approach of £(0,/3) to zero in 
the vicinity of T c can be studied by plotting InS versus ln(T c — T); we find that £ ~ (T c — T) x 
where the exponent x depends on A and lies between 0.4 and 0.6 for A between 1 and 2. This 
is consistent (up to subleading A-dependent corrections) with the value x = 1/2 characteristic 
of BCS theory. 

Figure 4 shows the scaled mass versus scaled momentum for several values of A, at a 
fixed value of (3. It is clear that £ decreases strongly as A increases from A = 1, suggesting 
that it may vanish for A larger than some critical A c . This possibility is examined in Figure 
5, which shows E(0,/?)/a versus A for various values of (3. We are not able to follow the S/a 
curves much below values of order 10 -5 due to numerical difficulties, but it seems reasonable 
to conclude that at a given temperature T, £ indeed vanishes for A > A c , where A c depends 
on T. For large (3 (low T), A c approaches a value just greater than 2. As the temperature is 
raised, A c decreases. At T = Appelquist et al. |2[ have found that in the limit A — > A c , the 
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zero-momentum mass vanishes according to 



£ 

— oc exp 

a 



-2vr 



VNJN - 1 



(11) 



We have explored the possibility of a similar behaviour in our model. As convergence near the 
critical point is very slow we have had to extrapolate from the calculated values of N to get the 
critical value N c for the corresponding temperature. Our results are shown in Figure 6 where 
we have plotted — ln(£(0, (3)/cx) vs. 1/ \JN C {T)/N — 1; we observe that for fixed temperature 
the curves approach straight lines as N approaches N c . This leads us to believe that indeed in 
this region the zero-momentum mass behaves like 



£(0,T) 

oc exp 



a 



-C(T) 



VN C (T)/N - 1 



(12) 



for some temperature-dependent function C(T) (which is, however, at the temperatures shown 
in Figure 6, considerably smaller in magnitude than the T = value of 2ir given in (|Tl|)). 

As mentioned in the previous paragraph, we find that as T — > N C {T) approaches a 
value just greater than 2. This is, of course, different from the N c found in Ref. Q, using an 
equation which ought to be the zero temperature limit of ours, and hence some further comment 
is required. In fact, while it is obviously true that the T — > limit of the full SD equation (||) 
must be the same as that of Ref. [0] , this is not the case after the instantaneous approximation 
has been made, leading to Eqns.(||) and Nevertheless, Eqn. (||) does reduce as T — > to an 
equation of similar form (for small £) to that in Ref. Q, but the numerical coefficient in front 
of the integral is a factor 1.5-2.0 times too small. Effectively this means that in comparing our 
results with those of Ref. [|| we should take our N as being roughly equivalent to the N of Ref. 
Q divided by this factor. This is the reason for the discrepancy in the N C (T = 0) values. 

These results enable us to obtain the phase diagram shown in Figure 7. There is a single 
critical line, such that for (N,T) below this line X ^ 0, and for (N,T) above it £ = 0. We have 
only shown the region N > 1, but it seems likely that the line approaches N = asymptotically 
as T — > oo. In this plot we have rescaled the critical line to match the zero-temperature results 
of Ref. §, namely N C (T = 0) = 3.2. 

As mentioned in the Introduction, the dimensionless ratio r = 2S(P = 0, T = 0)//cbT c is an 
important quantity, distinguishing between different mass-generation mechanisms. Our values 
of r are shown in Table 1 for N = 1, 1.5, 1.7, where for convenience we also list the corresponding 
values of S(P = 0, T = 0)/a and of kBT c /a. This table can be directly compared with Table 1 of 
Ref. |ll| which -it will be recalled- was obtained by solving a simplified equation in which the 
P-dependence of £ was neglected. The comparison shows that while our more exact equation (||) 
yields values of the mass £ which are about one order of magnitude smaller than those obtained 
in Ref. [11], the critical temperature T c is also correspondingly reduced, so that r remains with a 



value of order 10, in agreement with the value found in QllJ] , and also approximately independent 
of N. 

With an eye to the possible relevance to high-T c superconductivity, it is natural to wonder 
about the orders of magnitudes of the quantities appearing in our results, when expressed in 
physical units. We must emphasize, however, that relatively small changes in the kernel (which 
is after all only an approximation) can make rather large changes in £. It is clear that the 
scale of the model is set by the dimensionful quantity a, which has dimensions of (length)^ 1 or 
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(energy) in the system in which we have implicitly been working, namely that in which h = v = l 
with v the Fermi velocity in the condensed matter case (and that of light for QED3). If a 
is regarded as a freely adjustable parameter, then from Table 1 estimating ksT c ~ 10 _4 a for 
N = 2 (the value required by our high-T c model || ) , we would need to take a ~ 80 eV in order 
to obtain T c ~ W0°K as is required experimentally. It seems hard to understand how such a 
large energy could arise naturally. If we readjusted the "effective" N — T curve so as to agree 
with N C (T = 0) as found in Ref. we would obtain ksT c ~ lCP 3 a for N = 2, leading to a 
required a of order 8 eV . These different estimates merely underline, of course, the difficulty in 
making anything other than rather rough order of magnitude calculations as far as numerical 
values are concerned. A more reliable determination of E and T c would be obtained from finite 
temperature Monte-Carlo simulation of lattice QED3. 

A value of a in the region of a few eV is still much larger than typical Heisenberg exchange 
energies (recall that in our model || the gauge field arises in connection with the spin degrees 
of freedom in the original lattice Hamiltonian) . Nevertheless, it is possible to form an estimate 
of e 2 in terms of the parameters of the lattice model of ||, which shows that it is effectively 
enhanced, as follows. The lattice analogue of the fermion kinetic energy is the "hopping term" 
which enters with coefficient t. If the lattice fermion operators are rescaled by ta (where a is the 
lattice spacing) so as to get the correct dimensions of the fields in the continuum limit, and if 
space is then rescaled so as to obtain the (Dirac) kinetic energy with unit coefficient, the lattice 
U (1) coupling g becomes effectively replaced by gj (ta) 1 ^ 2 . An estimate of ta in such models may 
be obtained by noting that according to Baskaran et al. p4|] the maximum doping concentration 
n ma x ~ t/U, where U is the Hubbard repulsion. Since we may take U ~ a -1 (i.e. U — > 00 in 
the continuum limit) we find ta ~ n max , which has the empirical value of only a few percent. 
Assuming that the magnitude of g is set by the spin magnitude (1/2), and its length scale by 
the lattice spacing a, we obtain finally for the square of the effective coupling 

9 1 Hv , 
e ~ 7 (13) 

40. Tl max 

having reinstated Ti and v. The latter quantity can be conveniently found from the relation 
£ ~ Tiv/Yi for the correlation length £. Using £ ~ 30Aand S ~ 5ksT c we find v/c ~ 5 x 1CP 4 , 
which gives e 2 ~ few eV. Thus it is perhaps not impossible that such values could arise within 
the context of a model such as that of Ref. |J . 
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Tables 



N 


1 


1.5 


1.7 


£(P = 0,T = 0)/A 


1.93 x 1(T 2 


2.6 x 1(T 3 


9.0 x 10~ 4 


k B T c /A 


3.76 x 1(T 3 


5.2 x 10~ 4 


1.90 x 10" 4 


r = 2£(P = 0,T = 0)/fe B T c 


10.26 


10.0 


9.5 



Table 1: The zero-temperature and zero- momentum fermion mass and the critical temperature 
at a/ A = 1 with the ratio r for N = 1, 1.5, 1.7. 



Figure captions 



Figure 1: The photon propagator to leading order in 1/N 

Figure 2: Scaled dynamical mass E/a as a function of scaled momentum p/a for N = 1, and 
various scaled (inverse) temperatures /3a. 

Figure 3: Zero- momentum scaled mass versus scaled temperature for N = 1, 1.5. 

Figure 4: Scaled mass versus scaled momentum for /3a = 10 5 for N = 1, 1.2, 1.5. 

Figure 5: Zero-momentum scaled mass versus ./V at various (inverse) temperatures. 

Figure 6: Test for the behaviour given by Eqn. ( |i~2|) near the critical region. 

Figure 7: Phase diagram for spontaneous mass generation. 
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